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We obtained new periodic solutions in the problems of three and four point vortices moving on a plane. In the case 
of three vortices, the system is reduced to a Hamiltonian system with one degree of freedom, and it is integrable. In 
the case of four vortices, the order is reduced to two degrees of freedom, and the system is not integrable. We present 
relative and absolute choreographies of three and four vortices of the same intensity which are periodic motions of 
vortices in some rotating and fixed frame of reference, where all the vortices move along the same closed curve. Similar 
choreographies have been recently obtained by C. Moore, A. Chenciner, and C. Simo for the n-body problem in celestial 
mechanics [r>ll7lll7). Nevertheless, the choreographies that appear in vortex dynamics have a number of distinct features. 



Equations of motion and first integrals. Let us review briefly basic forms of equations and 
first integrals in the dynamics of point vortices on a plane (for a detailed discussion, see |U I14| , 
where, in addition, hydro dynamical assumptions required for the validity of these equations are spec- 
ified). 

For n point vortices with Cartesian coordinates (xi,yi) and intensities Tj, the equations of 
motion can be written in Hamiltonian form, 

where the Hamiltonian is 

n 

^ = --t^y iT M r i- r fii r i = (x i ,y i ). (2) 

i<j 

Here, the Poisson bracket is: 

U> 91 P, V PtT, fhu Fill: fir.- ) ' 1 ' 



Ti \dxi dm dyi dxi 

2 = 1 

Equations Q are invariant under the action of the Eucledian group E(2), so (in addition to the 
Hamiltonian) they also have three first integrals 

n n n 

Q = ^T iXi , P = Y,^iVu / = X^ r i(4 + yf)> (4) 

i=l i=l i=l 
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which, however, are not involutive, 

TV 

{Q,P} = J2^h {P,I} = -2Q, {Q,I} = 2P. (5) 

i=l 

Using these three integrals, it is possible to construct two involutive integrals, Q 2 + P 2 and I. 
Hence, according to the general theory jlUl I13j . we can reduce the number of degrees of freedom by 
two. Thus, in the particular case of three vortices, the system can be reduced to one degree of freedom 
and becomes integrable (Grobli, Kirchhoff, Poincare) |H 1141 Hoi I18| . while the problem of four vortices 
is reduced to a system with two degrees of freedom. Generally, the last problem is not integrable [20] . 

Effective reduction in the system of four vortices with intensities of the same sign was done by 
K. M. Khanin in P|. In that work two pairs of vortices were considered and for each pair action-angle 
variables were selected. The system of four vortices was then obtained as a perturbation of the two 
two-vortex systems. He proves (using the methods of KAM-theory) the existence of quasiperiodic 
solutions. As a small parameter, he takes a value inverse to the distance between two pairs of vortices. 

Reduction by one degree of freedom using the translational invariants P and Q was done by Lim 
in He introduced barycentric Jacobi coordinates (centered, in this case, at the center of vorticity), 
which have well-known analogs in the classical n-body problem in celestial mechanics [Sj- Note that 
even this (partial) reduction made it possible to apply some methods of KAM-theory to the problem 
of point vortices' motion [TT1 IT2] . 

When n vortices on a plane have equal intensities, one of the more formal methods of order 
reduction, which was described in [2] (see also USE]), proved to be most suitable. This method is 
based on using mutual variables representation of the equations of motion. As mutual variables, the 
squares of distances between the pairs of vortices and oriented areas of triangles were taken: 

Mij = (xi-xj) 2 + (yi-yj) 2 , A ijk = (rj - r k ) A (r k - r»). (6) 

Mutual commutation of such variables (which are introduced by E. Laura) leads to some Lie algebra. 
Here, the reduction procedure (more exactly, the last canonical step of it) is equivalent to a purely 
algebraic problem of introduction of symplectic coordinates on the orbits of corresponding Lie algebras. 
(For identical vortices the method of reduction based on the Fourier transformatio was introduced also 
in P.) 

Reduction for three and four vortices of equal intensity. Without losing in generality, 
we put Ti = Tj = l, P = Q = 0, then the moment integral / @ can be written as 

n 

J = kE M ^ ( 7 ) 

i<j 

where n is the number of vortices. Using the complex representation of the vortices' coordinates, 
z k = x k + Wki we obtain: 

n 

^ = IEV^^, (8) 

j=k 

where 6 k j is the angle between the vector from the j-th vortex to the k-th vortex and the positive 
direction of Orr-axis. 

The following propositions, defining the dynamics of a reduced system of vortices, can be obtained 
by direct computation: 

Proposition 1. For three vortices of equal intensity, evolution of the mutual distances © 
{assuming I = const) is described by a Hamiltonian system with one degree of freedom. In canonical 
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variables (g,G), the system looks as 



^ = HF' ^ = "^' ^ = -^ lnM i2^i 3 M 23 , (9) 



where M 12 = 4 U - g\ M 13 = 8G - I + 2^12 J U - G\ Gcosg, M 23 = 4 ( | - G 



- 2^12^ (| -G)g cos 5 - 

Proposition 2. For four vortices of equal intensity, evolution of the mutual distances is described 
by a Hamiltonian system with two degrees of freedom. In canonical variables (g,G,h, H), the system 
looks as 

9 = ~dG' G = ~~di' fc= W H = ~^h' (10) 
3V = In M 12 M 13 M 1A M 23 M 24 M U , 

where 



M12 = I-G + 2y/(I- H)(I -G) cash, M M = I - G - 2y/(I - H)(I - G) cos h, 
M 13 =I + G + 2^ (I - H)G cos(h + g), M 24 = I + G - H - 2y/(I - H)Gcos{h + g), 
M U = H + 2y/{H - G)G cos g, M 23 = H - 2y/(I - G)G cos g. 

Remark. These canonical variables are of natural geometrical origin related to representation of the 
equations of motion on a Lie algebra [21 E]. 

Absolute motion: quadratures and geometric interpretation. According to @, 
when Mij{t) are known, one needs to know the angles to determine the vortices' coordinates. 

It is obvious that only one of the angles is independent (in this case, we take 0\ 2 to be independent), 
the remaining angles are computed with the cosine theorem: 

/ M jk - Mij - M ik \ 

Oij + &ik = arccos 7==== — > i ^ J> k ^ L ( n ) 

\ 2yj2MijM ik J 

Evolution of Q\ 2 is obtained by quadrature ^1] : 

4tt#i 2 = M^ 1 + M^ 1 + M{^(6 - M^M 13 - M 1Z M^) (for three vortices), 
4tt# 12 = A^ 1 + Mf/ + M^ 1 + M2I+ (12) 
+ M± 2 l (8 - M{^M 23 - M^Mu - M 23 M 13 - M^Mu) (for four vortices). 

In the case of periodic solutions of the reduced system Q, (fT0|) . there is an interesting geometric 
interpretation of the absolute motion. 

Proposition 3. Let ^(t) be a periodic solution (of period T) of the reduced system, then 
1° there exists a frame of reference, uniformly rotating with some angular velocity Q, a about the 
center of vorticity in which each vortex moves along some closed curve £i(t); 

2° the rotational velocity Q a is given by (to within — -, p, qE Z): 



T 

b 



^a = h e l2 (t)dt. (13) 
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Fig. 1. The phase portrait of the reduced system in the case of three vortices. Relative choreographies 
corresponding to different phase trajectories on the portrait 

3° if the velocities Sl a and Sl Q = ^ are commensurable (i. e. ^ = ^, p,q 6 Z), then the vortices 

in the fixed frame also move along some closed curves; 

4° if any of the curves £i(i) can be superimposed by rotation about the center of vorticity by 
an angle, commensurable with 2ir, then there exists a (rotating) frame of reference such that the 
corresponding vortices move along the same curve. 

Proof. 

We expand the periodic function of period T in the right-hand side of (jl2|) in a convergent Fourier 
series 2 7rn 

4^i 2 = ^ a W e ! ~ t . (14) 



Integrating (|14j) while bearing (jllj) in mind, we come to a conclusion that the angles 6ij depend on t 
in the following way: 

e ij (t) = n ij t + g ij (t), (15) 

where Sly = + Tpj^T) Qij, Pij £ Z, while gij(t) = gij(t + T) are T-periodic functions of time. 
Substituting this into (jSJ), we see that locations of the vortices on a plane is given as follows: 

u k (t) = u k (t + Ti) G C, O = a (0) . (16) 



Zk (t) = Q±^ + Uk (ty n \ 



Hence, in £/ie frame of reference, rotating about the center of vorticity with angular velocity SI, all the 
vortices move along closed analytical curves, given by functions u k (t) £ C. 

The proof of 2°, 3°, 4° with is evident. ■ 

Remark. Proposition[2|is generalized without change to the case of arbitrarily many vortices, n, provided 
that j(t) is a periodic solution of the reduced system with 2n — 2 degrees of freedom. (For a more detailed 
discussion of reduction, see, for example, 
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Analytical choreographies. Now we show that the four- and three-vortex problems can have 
remarkable periodic solutions, when all the vortices follow each other along the same curve; such 
solutions are referred to as simple (or connected) choreographies. To underline the difference between 
choreographies in a fixed frame of reference and choreographies in a rotating one, these choreographies 
are called absolute and relative [7j. 

Theorem 1 If in the problem of three vortices of equal intensity, the constants of the 

integrals of motion, I and JF, satisfy the inequality 

-ln3 < ^je + ]nl < In 2, (17) 

then this motion is a simple relative choreography (see Fig. 
Proof. 

Since in the case of three vortices, the reduced system @ has one degree of freedom, all Mij(t) 
are periodic functions of the same period T. It can be easily shown that under the restriction (|17|). 
the orientation of the vortex triangle is unchanged, and there are times t± and t 2 such that M2z{t\) = 
= M 13 (t 2 ). (If JH| is not met, then Vt 3k such that M^k) < M ik (t) and My(k) < M jk (t).) 
Moreover, since with fixed H, I all are expressed through one of them (for example, M12), the 
following relations also hold true: 

Mi 3 (ti) = M 12 (t 2 ), M 12 (h) = M 23 (t 2 ). 

Since the equations are invariant under cyclic permutation of the vortices (as implied by equality of 
the intensities), we find that 

T 

h ~t 2 = g-rc, n G Z. 
Since the evolution equation for My is of the first order, we conclude that 

M l2 (t) = M 23 (*+!)= M 13 (t + ^ ) = f(t), 

or 



t + ^f) = M 13 (* + f ) 



where /(i) is some T-periodic function. 
Substituting this into (JSJ), we obtain 

^(t)=^(t + ^T^ e im , fc = 1,2,3, 

where n(t) is a T-periodic complex-valued function, determining the same curve, along which the 
vortices move, in the frame of reference rotating with angular velocity Vl. The phase portrait of the 
reduced system of the 3-vortex problem and the corresponding relative choreographies are given in 
Fig. [Q ■ 
The 4- vortex problem has an exceptional solution, given by quadratures — Goryachev's solution, 
where the vortices form a parallelogram at each instant of time [H]. As in the case of the 3-vortex 
problem, it is quite easy to show that 

Theorem 2 (|H). If in the J^-vortex problem, the vortices (of equal intensity) form a centrally 
symmetric configuration (a parallelogram), while the constants J4? and I satisfy the inequality 

-ln2 < ^J^ + lnl < -ln±|^, 

then the motion is a simple relative choreography. 

The corresponding relative choreographies are shown in Fig. |2 
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E = E r E= 1,22. = 4,09 E = 1 ,2552, fi = A A E = 1 ,35, S 7 = 5,945 
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o 



E = 1 ,364 w i' c , fi = 6,635 E = 1 ,37, Q = 6,88 £' = 1 ,6, S7 = 6,4 



Fig. 2. Relative choreographies in the 4-vortcx problem. 



Remark. The physical meaning of the inequality is as follows: when / is fixed, the type of the motion in 
the 3-vortex problem and in the case of Goryachev's solution is changed at energy values corresponding to the 
Thomson and collinear configurations. 

New periodic solution in the 4- vortex problem. Now we show that, aside from the chore- 
ography just described, the 4-vortex problem has at least one more choreography (different from 
Goryachev's solution). Consider the vicinity of Thomson's solution, i.e. the motion, where the vor- 
tices are located at the vertices of a square and rotate uniformly about the vorticity center ^j]- It is 
obvious that in the case of the reduced system with two degrees of freedom 1)10(1 . Thomson's solution 
is represented by a fixed point (more precisely, by six points corresponding to various arrangements of 
vortices at square's vertices). We consider one of the arrangenents (the rest are completely identical) 
with coordinates 

G = 0, H=\, fr = |. 
Let us find the normal form of the reduced system's (|1()|) Hamiltonian near this point. To do that, 




-0.25 

E= 1,215. n = 




0.25 -0.25 

5,11s 1 E= 1,22, (> = 5.205 




-0.25 o o:: 

F= 1,24. (> = 5,590 






-0.25 0.25 - 0.25 0.25 - 0.25 0.25 

E = 1,215, Q = 5,507 7 E= 1,26, ft = 6.0S* * E= 1,28, ft = 6.180 * 



Fig. 3. Relative choreographies corresponding to the new periodic solution in the three-body problem. 
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we start with a canonical change of variables: 

G=^±^, <? = arctgf, H=\ + Z-V% h = ^ + 8^y; 

Expanding the Hamiltonian in a series up to and including quadratic terms, we find 

^ = ln_2 + ^ 2 + J ^ (18) 

= h ( 3(x2 + x2) + 2 ^ (y2 + y2) ) ' 

where the expansion of starts with the third-order terms. Thus, the Hamiltonian ^tf 2 defines 
an integrable system with two incommensurable frequencies and has precisely two non-degenerate 
periodic solutions on each energy level 3P 2 = fi2 = const. These solutions are given by 

x = X = 0, y 2 + Y 2 = ^-h 2 , (19) 
y = Y = 0, X * + X 2 = ^fh 2 . (20) 




Fig. 4. Relative choreographies corresponding to the new periodic solution in the three-body problem in the 
frame of reference different from that used in Fig. [21 

According to the Lyapunov theorem [T^], these solutions are preserved under perturbations, hence, 
the complete system (|18[) in the vicinity of the fixed point also has two non-degenerate periodic 
solutions on each energy level. It is easy to show that one the solutions (corresponding to (|19|0 is 
identical to Goryachev's solution — during their motion the vortices are located at the vertices of 
a parallelogram. At the same time, the other solution (|2f)j) does not have such a simple geometric 
interpretation. 

Since equations (JJ) are invariant under a cyclic permutation of the vortices cr c (zi, z 2 , 23) = 
= (2:3, z\, z 2 ) and the eigenvalues of the Hamiltonian are different, it is easy to show that both 
periodic solutions are also invariant under a c . Thus, according to Proposition EJ in the appropriate 
frame of reference, all the vortices move along the same curve, i. e. both solutions correspond to 
simple relative choreographies. Figure El shows the relative choreographies that correspond to the new 
periodic solution of the reduced system (jlUj) . 
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Relative and absolute choreographies. Generally, for each periodic solution (of period T) 
of the reduced system Q, (|1U|) . it is possible to specify a countable set of rotating frames of reference, 
where the vortices move along closed curves. Indeed, the trajectories will remain closed in a frame 
rotating with velocity 

n' a = n a + ^, P ,qez (21) 

Nevertheless, the change (1211) with arbitrary p and q does not preserve the connectedness of the 
trajectories, i. e. in the general case, after moving to the frame of reference, rotating with velocity Q' a , 
a simple relative choreography decomposes into separate closed curves, along which the vortices move. 
To preserve the connectedness, the following criterion should be met. 

Proposition 4. Let a periodic solution (of period T) of the reduced system describe a connected 
relative choreography, the rotational velocity of the frame of reference being Q a , while the period of the 
motion of the vortices along the corresponding common curve being equal to mT. Then, if 

mp = knq, (22) 

where n is the number of vortices, and k £ Z is an arbitrary integer, the transformation (|21j) results 
in a connected choreography. 

Proof. 

For the solution in question, the absolute coordinates of the vortices can be presented in the form: 



k - 1 



iflat 



z k (t) = u[t+ ^j^mT ) e lUat , (23) 



where u(t) = u(t + mT) is a periodic complex- valued function (of period mT). 
Solving for Q a from (|2*T|) and substituting into (|2*3*|) . we get 

z k (t) =u(t+ ^mrj e -i # i* e in;* = Uk (t)e in '^. (24) 

If all the vortices move along the same curve, then their coordinates in the rotating (with velocity Q' a ) 
frame of reference are equal to Uk(t) and satisfy Uk+i(t) = Uk (t + ^jr-J, whence, taking into 



n j 

account, we obtain (|22|). 



The relation (|22|) is a sufficient but not necessary condition for the choreography to be connected. 
If the curve along which the vortices move has additional symmetries, then, besides the p, q that meet 
the condition ()22|) , there are more velocities of the form (|21(l , which result in connected choreographies 
(see below, for Goryachev's solution). 

It is interesting to note that using the transformation (|21j) . some choreographies can be "disen- 
tangled" — for example, Fig. |I] shows choreographies, corresponding to the periodic solution ()2(Jj) . in 
the frame of reference, different from that used in Fig. 13 their rotational velocities differ by 

n'-n = l 2 f. 

As it was shown above (see Proposition EJ), if, for a relative choreography, the period T of the 
reduced system's solution is commensurable with the rotational period T a = j^- of the frame of 

* 'a 

reference, then in a fixed frame of reference, all the vortices move along closed (and, usually, different) 
curves. 

Let us consider in greater detail the existence of absolute choreographies in the three- and four- 
vortex problems. According to what was said above, any relative choreography, corresponding to the 
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Fig. 5. The rotational angular velocities, fl\ and fl 4 , of relative choreographies are presented as functions 
of energy. The points where the graphs intersect the Ox axis correspond to the absolute choreographies of the 
three vortices shown in the bottom (Figs, b and c). The heavy lines denote the basic angular velocity fij , 
corresponding to the simplest relative choreography (shown in Fig. a) , and the frequency of the periodic solution 
of the reduced system 101, Q, = —. 




Fig. 6. The rotational angular velocities, fij an d of relative choreographies are presented as functions 

of energy. The points where the graphs intersect the Ox axis correspond to the absolute choreographies of the 
four vortices shown in the bottom (Figs, b and c) . The heavy lines denote the basic angular velocity > 
corresponding to the simplest relative choreography (shown in Fig. a), and the frequency of the reduced sys- 
tem JUJ), n = 
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periodic solution (of period T) of the reduced system (see @, (|10|0 . closes in time mT, m € N. During 

(k) 

this time interval, the vortices pass one and the same relative configuration m times. Let Sim be the 
rotational velocities of the frames of reference related to these choreographies. 

As it was shown above, all the solutions of the reduced three-vortex system, for a fixed D 
and Ef < E < Eq (where Et and Eq are the energies corresponding to the Thomson and collinear 
configurations), describe connected relative choreographies. Moreover, there is a frame of reference, 
where the choreography closes in (the shortest possible) time T (see Fig.[5J); the corresponding angular 
velocity is denoted as , its graph is shown in Fig. [5J The angular velocities of the other relative 
connected choreographies are now given by 

Q$(E) = nf\E) + ||fio(-E), where m£N,kZ; (25) 

here 3k and m are coprime numbers, and Qq(E) = ^ . The velocity Sim corresponds to the 

T(E) 

choreography that closes in time mT. 

Obviously, the absolute choreographies are defined by the solutions of the equation 

n$(E) = 0, (26) 

where k, m are fixed, while E is unknown. Figure El shows the graphs of velocities Qm\E) together 
with the solutions of equation (|26|). as well as the corresponding absolute choreographies (the simplest 
choreographies in a three- vortex system). Generally, there is a countable set of absolute choreographies 
with different m, k. 

Indeed, consider the function 

f a (E) = nf ] (E)+an (E), (i£t, 

for which the equality f a (Ec) = wc < holds (see Fig. [SJ. Since £lo(E) > 0, there exists a number a* 
such that when a > a*, the function f a (E) has at least one zero in the interval [Et, Eq\. It is clear 

that the interval [a*, +oo) contains an infinite number of rationals of the form a = — , where 3k 
and m are coprimes. The dotted line in Fig. |S] shows the curve f a *(E) and the relative choreography's 
rotational angular velocity £l a , calculated from (|13|). This choreography is the simplest disconnected 
choreography and is remarkable for the fact that the values at the ends of the interval, Q a (Er) and 
0, a (Ec), are equal to the rotational angular velocities of the Thomson and collinear configurations. 

For Goryachev's solution in the four-vortex system, the reasoning is similar but slightly modified. 
First of all, one can show that the simplest connected choreography closes in time 2T, while the 
vortices in this case pass one and the same relative configuration twice (i. e. velocities correspond 
to disconnected choreographies). The graph of one of the corresponding angular velocities, which 
we denote as is given in Fig. ??. In this case, due to the symmetry of the curve related to 

the choreography the frequencies corresponding to connected choreographies must satisfy the 

relation different from (I25|) . 

n%l(E) = l4 0) (£) + m Q o(E), where m is odd, kZ; (27) 

here k and m are coprimes, and Q,q{E) = ^ 7r , where T is the period of Goryachev's solution to the 

T(E) 

reduced system H10|) . This choreography closes in time 2mT. 
As above, the equation 

42(^) = o 

defines the absolute choreographies. Similar to the three-vortex problem, one can show that there is 
a countable set of absolute choreographies, described by Goryachev's solution with different m and k. 
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Stability. As a conclusion, let us discuss the stability of the specified periodic solutions. 

Since the three- vortex problem is integrable, all the solutions of the reduced system Q are 
periodic and stable. Yet it is easy to show that any (absolute or relative) choreography in this 
problem is neutrally stable under perturbations of the vortices' positions in the absolute space. 

In the four-vortex problem, due to its non-integrability, the relative choreographies can be (expo- 
nentially) unstable. Yet, if a periodic solution of the reduced system ifTUjl is stable, the corresponding 
choreographies are also neutrally stable in the absolute space. The numerical analysis of the multi- 
pliers of the periodic solutions that describe choreographies presented in Fig. ?? b, c shows that they 
are (exponentially) unstable. 

We gratefully acknowledge support of the Russian Foundation for Basic Research (grant 
JNffi 04-05-64367), the Foundation for Leading Scientific Schools (grant N°- 136.2003.1) and CRDF 
(grant RU-M1-2583-MO-04). 
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